Semiclassical quantization of skipping orbits 
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We propose a simple description of the spectrum of edge states in the quantum Hall regime, in 
terms of semiclassical quantization of skipping orbits along hard wall boundaries, A = 2TT{n + ^)£%, 
where A is the area enclosed between a skipping orbit and the wall and £b is the magnetic length. 
Remarkably, this description provides an excellent quantitative agreement with the exact spectrum. 
We discuss the value of 7 when the skipping orbits touch one or two edges, and its variation when 
the orbits graze the edges and the semiclassical quantization has to be corrected by diffraction 
effects. The value of 7 evolves continuously from 1/2 to 3/4. We calculate the energy dependence 
of the drift velocity along the different Landau levels. We compare the structure of the semiclassical 
cyclotron orbits, their position with respect to the edge, to the wave function of the corresponding 
eigenstates. 

PACS numbers: 



I. INTRODUCTION 

The edge states play a crucial role for understanding 
the integer and fractional quantum Hall effects. Their 
description has been introduced in the seminal paper by 
Halperin.— This picture has then been elaborated by 
Buttiker— and a review can be found in ref. It ap- 
pears often convenient to picture qualitatively these edge 
states in terms of skipping cyclotron orbits. But the link 
between the full quantum mechanical treatment of the 
states and this qualitative picture is missing (see however 
ref. 0). Here we propose an extensive development of 
this picture and show how the semiclassical quantization 
of these orbits leads to a qualitative and even quantita- 
tive description of the edge states energy levels. 

We consider a free electron (mass m, charge — e) mov- 
ing in a ribbon infinite along the y direction and bounded 
along the x direction. A magnetic field B is applied along 
z. As in ref. we consider the situation where the 
confining potential consists in an abrupt potential well 
of infinite height. This is known not to be the correct 
situation in the two-dimensional electron gas of GaAs 
heterostructures, where the confining potential is rather 
smooth at the scale of the magnetic length £b- How- 
ever, we believe that the case of the abrupt potential is 
interesting in itself and may be relevant to other related 
situations exhibiting edge states. For example, graphene 
ribbons have sharp boundaries which must be modeled 
with sharp potentials.^'^ We comment the case of smooth 
boundaries at the end of the paper. Far from the edges, 
the energy levels are given by i?„ = (n + l/2)hLUc where 
Wc = eB/m is the cyclotron frequency. This Landau 
quantization can be obtained quite easily from the Bohr- 
Sommerfeld quantization rule that we recall below. The 
goal of this work is to describe the semiclassical motion 
of the electron near one edge, described here by a "hard 
wall" , that is an infinite potential well. We consider the 
vicinity of the edge located at cc = 0, assuming first that 
the second edge is far away {x — > —00). 

The problem is solved semiclassically by quantization 



of the action. In the appropriate gauge, the well-known 
Landau problem is related to the problem of a one- 
dimensional oscillator. In the presence of a sharp edge, 
the problem to be solved is the one of an harmonic os- 
cillator in the presence of an infinite potential well. This 
problem is solved by quantization of the semiclassical ac- 
tion S{E) — 2TT{n + 'y)h, where 7 is related to a so-called 
Maslov indexj^ For the free oscillator, 7 = 1/2 corre- 
sponds to the sum of two contributions % — 1/4 of the 
two turning points. In the presence of the potential well, 
7 evolves from 1/2 to 3/4 when the guiding center Xc of 
the cyclotron orbit (the center of the harmonic oscillator) 
approaches the wall. A form of the continuous variation 
Jn{xc) for a given n has been recently obtained.^ 

In this paper, we give a very simple description of the 
edge states spectrum in terms of quantization of skipping 
orbits. This image, currently used in the literature or in 
pedagogical presentations, has curiously never been de- 
scribed in details (see however ref. jlj). Yet, it leads to 
a number of results which to our knowledge have never 
been discussed. In the next section, we recall the map- 
ping, in the Landau gauge, to a one-dimensional problem 
of a harmonic oscillator and we calculate the action S{E) 
of this oscillator. In section IIIIl we give a complete pic- 
ture of the evolution of the energy levels in terms of the 
quantization of the area of skipping orbits. The well- 
known quantization of closed orbits can be extended to 
the case of skipping orbits. Then, their area depends on 
the distance Xc to the wall and must be quantized as 



A{R,Xc) = 27r(n + 7)^| 



(1) 



where £b is the magnetic length, R is the cyclotron radius 
and Xc is the position of the guiding center with respect 
to the wall. This well-known quantization rule for closed 
orbits appears to be also valid for open but periodic skip- 
ping orbits. Using the same method, we calculate in sec- 
tion IIVI the full spectrum in the case a ribbon, when the 
magnetic length is of the order of the width of the ribbon 
so that the two edges have to be considered. Then we 
conclude in section |Vl with a comparison with the case 
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of a smooth potential. 



II. MAPPING TO A ONE-DIMENSIONAL 
OSCILLATOR 

The problem to be solved is described by the Hamilto- 
nian 

where the potential V{x) describes the edge of the sample 
along the y direction. We choose V{x) =0 when a; < 
and V{x) = oo when a; > 0. Using the Landau gauge 
A = {0,Bx,0), the corresponding eigenvalue problem 
reads (— oo < x < 0) : 

[-|^TT + ^^^c(^c - xfmx) = E^ix), (3) 
2m ax'' 2 

with the constraint V'(O) = 0. The center Xc of the os- 
cillator, is related to the ky component of the wave vec- 
tor which is a good quantum number: Xc ~ kyi'^ where 
£b = \ffifeB is the magnetic length. The action S along 
a closed trajectory is given by 5 = f Pxdx — 2 J p^dx 
where, far from the edge (region I in Fig. [T]^ 

S = 2 / y^2mE — w?uj'^{x — Xc)'^dx 

2'kE o 
= = irmujcR . (4) 

R{E) = y/2E/mujl is the cyclotron radius of the classical 
trajectory and Xc ^ R are the positions of the turning 
points. The semiclassical quantization of the action 



Introducing the angle such that cos6' = Xc/R, the ac- 
tion can be rewritten as 

E 1 
S{E,x^) = — [261 - sin 26*] ^ -mtJc-R^[26' - sin 261] . (8) 
cjc 2 

We give in the next section a very simple interpretation 
of this angle 0. The energy levels En{xc) are still given by 
quantization of the action ([5]) which now depends on the 
position with respect to the wall. When Xc > —R, the 
factor 7 is equal to 3/4, because it results from different 
matching conditions at the two turning points. At the 
left turning point ji — 1/4, while at the right turning 
point, the vanishing of the wavefuntion implies 7^ = 1/2, 
so that 7 = 7; -I- 7j. = 3/4. 7 evolves between 1/2 and 
3/4 when Xc — —R. 
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FIG. 1: ID harmonic oscillator centered on Xc, with an infi- 
nite potential well inx — 0. a) The guiding center is inside the 
sample, b) The guiding center is outside the sample. Three 
distinct regions have to be considered : I : the energy levels are 
not affected by the edge and 7 = 1/2, //; the trajectories hit 
the edge and 7 = 3/4. In between, the phase factor 7 evolves 
between 1/2 and 3/4. It is 2/3 when the cyclotron orbit just 
grazes the wall (III). 



S{E) ^2TT{n + j)h (5) 

leads to the energy quantization, 

(n + -f)hu!c (6) 

where the value 7 = 1/2 is not given by the semiclassical 
quantization rule and results from the matching of the 
wave function at the two turning points. 



III. QUANTIZATION OF SKIPPING ORBITS 

A. Quantization of the area 

The quasiclassical Bohr-Sommerfeld quantization rule 
for skipping orbits has been discussed by Beenakker and 
Van Houten^. Since the motion along the x axis is peri- 
odic, this quantization rule can be written as 



When the cyclotron orbit approaches the edge, that is 
when the center Xc of the cyclotron orbit becomes larger 
than —R (energy regions Ila and lib in Fig. [ij, the 
turning points are located at xi = Xc — R and X2 = 0. 
The action now explicitly depends on Xc and it is given 
byl2 



S{E,Xc) = 2 / ^2mE - m?ujl{x - x^)'^dx 

= — [f-arcsin^-^Vl-^M-(7) 

LJc ^ ti IX 



S ^ j) ■ dx ^ 2TT{n + j)h . (9) 

where p = mv — eA. The trajectories are now open and 
the notation f means that the integral is taken along one 
period of the motion. The gauge for the vector potential 
A must be chosen such that is periodic. The sim- 
plest choice is A^ = 0, so that px = mvx- The classical 
equation of motion for the x component of the velocity 
is Vx = —i^c{y — yo) where yo is an arbitrary position. 
Therefore the quantization condition ^ becomes 
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(a) 



(b) 



FIG. 2: Semiclassical skipping orbits a) when the guiding cen- 
ter is inside the sample (region Ila in Fig. Qp; b) when it is 
outside (region lib in Fig. QP. The dashed area A is quantized 
as A = 2Tr{n + "/)£%. 



dx — —eB f{y~ ya)dx ~ 27r(n + "f)h . (10) 



The integral is nothing but the area A enclosed between 
one arc of the periodic orbit and the wall (see Fig. [2|). 
Therefore we can generalize the familiar quantization rule 
([T]) of the area A to the case of skipping orbits. This area 
can be parametrized by the angle 6 shown in Fig. [2] and 
defined hy Xc = R cos 9. We have 



R 

AiE,Xc) = — [261 - sin 261] 



(11) 



which is precisely the same equation ([8|) as obtained in 
the ID picture. Then the quantization of the area A 
reads 



A{R,x,)^2Tr{n + -f)el 



(12) 



so that the angle can be used to parametrize the solu- 
tions {0 = tt: the orbit just grazes the edge, Xc = —R- 
9 = 7r/2: the guiding center of the orbit in precisely on 
the edge. 9 < tt/2: the center of the orbit stands outside 
the sample). From equations (|llll2p . we obtain 



r:' 



47r(n + j)£' 



29 - sin 26* 
and the energy levels are given by 

i?2 



E — fkJr 



2l\ 



huJc{n + 7) 



2tt 



29 - sin 26* 



(13) 



(14) 



The cyclotron radius R is related to the position Xc 
the guiding center: 



of 



= i?cos = 



47r(n + 7) 
26* - sin 26* 



cos 



(15) 



so that the dependence En{xc) is simply parametrized 
by the angle 0. However, the main complexity of the 
problem comes from the fact that 7 is not a constant. 
It is fixed to the value 1/2 far from the edge, but on 
the other hand, for skipping orbits, it reaches the value 
3/4. Therefore, from the quantization condition ([T^. we 
obtain two branches (Fig. [3]). 



B. Spectrum 




FIG. 3: Energy levels obtained from semiclassical quantization 
of the area A{E,Xc) — 2-K{n+'))Pg. Inside the sample (xc « 
—R), 7 = 1/2. For skipping orbits (xc » ~R), 7 = 3/4. 
The dashed line corresponds to Xc — —R, the case where the 
cyclotron orbit just touches the edge. In this case, we have 
found that 7 = 2/3 and the position of the energy levels is 
marked with small dots along the dashed line. The energy 
levels are plotted in units of hwc and the distance Xc is plotted 
in units of £3. 



In the intermediate region, when the cyclotron orbit is 
very close to the wall, that is when Xc — —R, 7 varies 
continuously between 1/2 and 3/4. This regime has been 
studied recently within a WKB approach.^ In particular, 
when the cyclotron orbit strictly touches the wall Xc = 
—R, it has been found that the parameter 7 = 2/3. This 
factor comes from a contribution 1 /4 on the left side and 
a very peculiar and new contribution 5/12 from the right 
side, which, to our knowledge has never been studied, at 
least in this context. Moreover, in ref. Q, we have found 
an interpolation formula for -fn{xc) for a given value of 
n. It is given by 



AX 



1 l + 3e 
^^(^-) - 2 l + 2e'4x 



(16) 



where X ={2n + 4/3)'^/^{Xc/iB + ^2n + 'i/3) and A ~ 
3.5. This expression can be extended by transforming it 



4 



into a function of energy and Xc : 7(-E, a^c) is still given 
by Eq. (HH), with X = {2E /hwcY/^{xc/ Ib + y/2EjhJl). 
It can be actually decomposed in the form 

l{E,x,) = ]^+-fr{E,x,) (17) 

since it is known to be the contribution of two terms 
corresponding respectively to the left and to the right 
turning points. We have 



7^(^'^c) = i (18) 

Given these expressions and the implicit equation (1121) , 
the full spectrum is obtained in Fig. |4l 

The scenario when the cyclotron orbits approaches the 
edge is the following. When Xc ^ —R, that is far from 
the edge, the cyclotron radius is i? = iB^/^n + 1. When 
the distance between the orbit and the edge becomes of 
order of the magnetic length the energy and the 
cyclotron radius start to increase t o reach the values 
En = {n + 2/3)fia;c and R = £B\/2n + 4/3 when the 
cyclotron orbit just touches the edge. Then the energy 
and the cyclotron radius continue to increase as shown 
in Fig. m 
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FIG. 4: Black curves: Energy levels obtained from semiclassi- 
cal quantization of the area A{E,Xc) = 2-K{n + 'yn{xc))t%, 
where we have used the expression 116\) for 7n{a;c). The 
dashed line corresponds to Xc = —R, the case where the cy- 
clotron orbit just touches the edge. On this line, 7 = 2/3. The 
dotted lines are the approximation of a constant j — 1/2 or 
3/4, see Fig. [5l The large dots indicate special points where 
the wave function can be easily obtained from the solution of 
the harmonic oscillator in free space (see Figs. |g|_7|). When 
Xc — —R, 7 varies continuously between 1/2 and 3/4 so that 
the spectrum is continuous. The energy levels are plotted in 
units of hwc and the distance Xc is plotted in units of £b . 

It is also interesting to introduce the position h of the 
extremum of the cyclotron orbit (Fig. [2]), that is the 
position of the left turning point in the ID picture. It is 
h — —2R when the cyclotron orbit just touches the edge 
and it varies to /i — > when Xc — > 00 . A simple geometric 



picture shows that h = Xc — R = R{cos6 — 1) < 0, that 
is, using [TH 

h = Ricose-i)^iB^[^^^ico.e-i) (19) 

In Fig. [5l we plot the energy as a function of the position 
h. Of course, Xc can increase to infinity and h stays 
confined to the inside of the sample {h <0). (cf. inflexion 
point). 
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FIG. 5: Energy levels En(h) obtained from semiclassical quan- 
tization, h is the position of the extremum of the cyclotron 
orbit. The dashed line corresponds to h — ~2R, the case 
where the cyclotron orbit just touches the edge. On this line, 
7 = 2/3. The dotted line corresponds to h — —R, that is 
Xc = 0. The large dots indicate special points where the wave 
function can be easily obtained from the solution of the har- 
monic oscillator in free space (see figures [gj^j. The energy 
levels are plotted in units of hiOc and the distance h is plotted 
in units of I b 

The large dots marked in Figs. I4I5I correspond to sim- 
ple cases where the energy and the wave function are 
easily known. For these special points, where the en- 
ergy is the same as in free space (3/2,5/2,7/2 x hwc), 
the wave function is also the same as in free space but 
must vanish in a; = 0. The wave functions in free space 
are well known to be related to the Hermite functions. 
Therefore the edge must coincide with a zero of these 
Hermite functions. For example the points B and H cor- 
respond to antisymmetric wave functions, that is to en- 
ergies En = E2P+1 = huJc{2p -t- 3/2). In Fig. [SI we have 
shown the evolution of the normalized squared wavefunc- 
tion I V'(a;) p for increasing values of Xc. In Fig. [71 we show 
three wave functions with energy 5/2?uJc, having respec- 
tively two, one and no zeroes. In these two figures, one 
sees that the extension of the wave function is given by 
the extremum of the classical skipping orbit. 

C. Drift velocity 

We now calculate semiclassically the drift velocity 
along an edge state at a given energy E. The cyclotron 
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FIG. 7: Evolution of the state of energy 5/2fia;c when Xa in- 
creases, corresponding to the points (E), (F) and (C) shown 
in Figs. [i[5l The extension of the wave function is given 
by the extremum of the semiclassical orbit. Distances are in 
units of Ib. 



clotron orbit is given hy ve — ^^cR — \/2£'/m. The 
length of the skipping orbit being 2R9 (Fig. [2]), the pe- 
riod T is given by 



FIG. 6; Evolution of the ground state wave function (n — 0) 
when the position Xc of the guiding center approaches and 
crosses the edge. The four cases correspond to the points (A), 
(B), (C) and (D) shown m Figs. [i|5l The position of the 
turning point (at distance \h\ from the edge) is marked with a 
dot. Is is seen that the extension of the wave function is given 
by the extremum of the classical skipping orbit. The distance 
x IS plotted in units of Ib ■ 



T = 2R6/v = — 

UJc 



(20) 



Classicahy, the drift velocity v'jl along the y direction can 
be easily obtained from simple geometry (Fig. [5]), the 
distance between two successive hits on the edge being 
2i?sin6l: 



radius is i? = ^j2W[mj2^ and the velocity along the cy- 



2i?sin6' 
T 



sin Q 



VE- 



(21) 
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It is interesting to compare this value to the drift velocity 
obtained from the energy 



Vd 



dE _ I dE 
hdky eB dxc 



(22) 



which is the correct result, beyond semiclassical approxi- 
mation. The derivative can be calculated from Eqs. ([HI 
[T5|) . and one recovers the classical expression (PT|) pro- 
vided 7 is a constant. However in the region where the 
cyclotron orbit is near the edge Xc — — -R, 7 varies con- 
tinuously between 1/2 and 3/4. A better evaluation of 
the drift velocity is obtained in the WKB approxima- 
tion which accounts for the variation of 7 (Eq. [T6l) . The 
energy levels are given by the two equations 



47r[n-|-7„(xc)] 
20 ~ sin 26 



i?cose' (23) 



Since E = hujcR^ /2(.\, the drift velocity is 



WKB 



— —-LjR— (24) 
eB dxr dxr 



By differentiating Eqs. 
obtain 



WKB 



i^cR- 



2lTUJ, 



, and eliminating dR/d9, we 



sm 



c'^B' 



(1 — cos 26) dxc 



(25) 



The derivative is non zero only in the vicinity Xc — 
— i?, in which case the angle 6 is very close to tt. There- 
fore, in a very good approximation, we obtain: 



WKB 



n2 din 

'^^dx. 



(26) 



The dependence v^^^ixc) is plotted in Fig. ^ and 
is compared to the numerical fully quantum calculation. 
The drift velocity, which is zero inside the sample, starts 
to increase when the cyclotron orbit touches the edge. 
The WKB approximation is excellent except close to the 
point where the classical cyclotron orbit just touches the 
boundary. When the skipping orbit gets closer and closer 
to the edge, the energy increases and the classical approx- 
imation (PT|) becomes excellent. 

Fig. ini represents the drift velocity normalized to the 
Fermi velocity, as a function of the energy along a given 
edge state. The drift velocity ultimately saturates to- 
wards the Fermi velocity at high energy. 



IV. TWO EDGES 

The description of a ribbon with two edges is straight- 
forward when the two edges are sufficiently far apart com- 
pared to the cyclotron length. Here we consider the situ- 
ation where this is not necessarily the case, that is when 
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FIG. 8: Variation of the drift velocity Vd (in units of h/mlB ) 
with the position Xc for the three lowest energy levels. Red 
dots: exact numerical calculation; full line: result S26\) of the 
WKB calculation; dashed line: classical drift velocity \21]) . Xc 
is in units of £b- 
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FIG. 9: Variation of the drift velocity with the energy along 
different Landau levels. The drift velocity is normalized to the 
Fermi velocity ve corresponding to energy E. The full lines 
are the WKB calculations and the dots are the exact results. 



the distance d > Q between the two edges is of the or- 
der of a few magnetic lengths £3- The spectrum with 
two edges obtained numerically for d = Sis is shown in 
Fig. (fTO| . and clearly exhibits three different regions. 
We now give a full semiclassical description of this spec- 
trum, considering these three different cases. Regions I 
and II have already been discussed and correspond either 
to a free cyclotron orbit or to an orbit skipping along one 
boundary. The new interesting case is the region HI for 
which a cyclotron orbit touches the two boundaries. 

We first show that in this case the area A to be quan- 
tized is the area of a circular orbit cut by the two bound- 
aries (area HI in Fig. [T^ . This may not seem a priori 
obvious since this area is not bounded by a classical tra- 
jectory. Actually a periodic trajectory, the arc irf in 



7 




FIG. 10: Landau levels spectrum calculated numerically for a 
ribbon of finite width d — 8£b- The two vertical lines indicate 
the position of the edges and the two parabolas indicate the 
positions Xc for which the classical orbits touch the edges, Xc ~ 
~R and Xc = — d + R. Xc is written in units of Ib- 




A" 



FIG. 11: In the presence of two walls, the area delimited by a 
periodic trajectory in the area J\! , but the area A to be quan- 
tized IS A = A + A" . 



Fig. [TTJ encloses an area A' smaller than A. Let us re- 
turn to the argument developed in section (III. A). The 
Bohr-Sommerfeld quantization rule states that the inte- 
gral of the velocity along the trajectory irf has to be 
quantized: 



Vxdx — 27r(n + 



(27) 



i and / are respectively the initial and final points of the 
periodic trajectory and r is the point where the trajec- 
tory bounces on the second boundary (Fig. [TT|) . Now 
the velocity must be calculated with caution. Along the 

trajectory ir , it is given by Vx — —(^c{y — Vo)- Then 
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FIG. 12: Particular trajectories showing the four different 
possibilities described in Fig. S10{) . The spectrum is obtained 
from quantization of the shaded areas. 



after the bouncing along the second wall, it is now given 
by Vx = —ujc{y + Urr' — Vo) whcrc i/rr' is the distance be- 
tween the bouncing point r and the point r' which is the 
next intersection between the fictitious cyclotron orbit 
and the second boundary. Therefore we have: 



1 



V^rdx 



{y - yo)dx - 



Hrr'dx (28) 



The first integral on the right side is the area delimited 
by one period of the motion [A' in Fig. [TT|) and the sec- 
ond integral is the shaded area {A" in Fig. [TTT) . The sum 
of these two areas A — A' + A" is indeed the total area 
delimited by the free cyclotron orbits and the boundaries 
(Fig. [UlII). 

Defining d > as the distance between the edges, this 
area A is now given by 

A{E, Xc) = — [29r - sin 29r - 201 + sin 261/] (29) 

where 9r and 9i define the position of the cyclotron or- 
bits with respect to the two edges (see Fig. [121 We 
have cos 9r = x^jR and cos0; = (d -f x^jR, where 
—d < Xc < ^ when the guiding center is inside the 
sample. The energy levels are semiclassically given by 
the quantization ((T2|) of the area A{E, Xc) given by (|29l) . 
where the index 7 depends on the geometry of the orbit. 
It is the sum of two terms, 7 = 7/-!- 7r, where 7/,^ = 1/4 
in free space, 7/,^ = 1/2 for a skipping orbit, 7;,^ = 5/12 
when the cyclotron orbit just touches a wall. Between 
these different values, 7 varies continuously. In the gen- 
eral case, we have obtained the value of 'y{E,Xc) from 
its decomposition explained above (Eq. [T7)) . Its value is 



given by 



20 F 



7 = 7/ + 7r 



with 



Ir.l 



1 1 + 46-4-^- 

4 r 



2e 



(30) 



2E 



1/6 



— ) (a;, + ^2E/nuj,) 



2E 

husc 



1/6 



The function "f{E,Xc) is shown in Fig. [T3]as a function 
of the energy and the position Xc in the ribbon. Note 
that in the hmit where the ribbon is narrow d << i?, 
that is in the high energy regime III, we recover straight- 
forwardly that the area is now A{E, Xc) — 2dR, so that 
the quantization of this area gives R = Tr{n + ^)£'^/d and 



En — 



o ^T-, with n' — 
ing to the two reflections on the edges 
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FIG. 13: Bottom: dependence "/{E, Xc) versus Xc for various 
fixed energies (E/huJc = 2,4,6,8,10,12, from bottom to top 
curves). We have indicated the special values |i |i f ! f i 1 1 
corresponding to the different regions shown on the upper di- 
agram. 

The spectrum obtained from semiclassical quantization 
of the area P5|) with fixed values of 7 = 1/2,3/4, 1 cor- 
responding to an free orbit, an orbit touching one or two 
edges, is displayed in Fig. 1141 The approximation is 
quite good but there are discontinuities corresponding to 
Xc = —R and Xc — —d + R. In Fig. [T5l the spectrum 
is obtained from quantization of the area, with the ap- 
propriate value of 7 obtained above (Eqs. I30l3ip . We 
obtain a perfect quantitative description of the full nu- 
merical spectrum. 
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FIG. 14: Full curves: Spectrum obtained by semiclassi- 
cal quantization of the area h29{l . with fixed values of 7 = 
1/2,3/4,1 corresponding respectively to free cyclotron orbits, 
orbits touching one or two edges. Dotted curves: Exact spec- 
trum obtained by numerical calculation. The semiclassical ap- 
proximation is quite good except when the classical cyclotron 
orbit approaches the edges. The two parabolas correspond to 
Xc ~ —R and Xc ~ — d + R. 



V. CONCLUSION 

We have provided a semiclassical treatment for the po- 
sition dependence of the edge states energy levels in the 
presence of an abrupt infinite potential. This full spec- 
trum may be obtained from the Bohr-Sommerfeld quan- 
tization of the area of cyclotron orbits. The orbits do not 
need to be closed, and the quantization is obtained in all 
cases, where the orbits hit one or two edges. We pro- 
vide a simple expression for the mismatch factor 7 valid 
for all energies and positions with respect to the bound- 
aries. The situation of an abrupt potential corresponds 
to a physical limit where the range of variation of the 
potential at the edge is much smaller than the magnetic 
length £b. In the case of a smooth potential, the correct 
description corresponds to the adiabatic approximation 
where the energy levels simply follow the potential V{xc) 
at the edge: En{xc) = {n -\- l/2)hu}c + V{xc), as shown 
in Fig. [ini An important difference with the case of the 
abrupt potential is that here the energy profile is exactly 
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pends on n, it starts to increase at a distance of order 
\/2n + \£b from the boundary (Compare Figs. |4] and 
[T6| . Another important difference is that, for an abrupt 
potential, the maximal value reached by the drift veloc- 
ity is of order of the Fermi velocity ve- If one considers 
a soft potential of the form muj'^^x + a;o)^/2, where it is 
usually assumed that lo < ljc (this corresponds to the 
approximation IbV < fvjjc), the maximal velocity is of 
order u^joj^/cj^, much smaller than ve- 

In conclusion, we have shown that the semiclassical 
picture of quantized skipping orbits leads to a quantita- 
tive description of the edge states energy spectrum. We 
believe that this quite simple description, not only has a 
pedagogical interest, but may allow the study of physi- 
cal quantities not very much discussed in the literature, 
like the drift velocity. We believe also that it can help 
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FIG. 15: Full curves: Spectrum obtained by semiclassical 
quantization of the area f29i) . with 7 given by Eqs. \3U\3Vil . 
Dotted curves: Exact spectrum obtained by numerical calcu- 
lation. This semiclassical WKB approximation is now perfect 
even when the classical cyclotron orbit approaches the edges. 
The two parabolas correspond to Xc = —R and Xc = — d -|- R. 



the same for all levels. In particular, the drift velocity is 
independent on n and is simply given by 



1 dV 

eB dxc 



(32) 



and it starts to increase when the potential increases, 
while for the abrupt potential, the drift velocity de- 



FIG. 16: Energy levels in the case of a smooth confining po- 
tential V{x). In this case, the energy levels following simply 
the potential profile: En{xc) ~ {n-{-l/2)hujc-\~V{xc), compare 
with Fig. 



for the description of more sophisticated problem like the 
structure of edge states in graphenei^ 
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